Preprint typeset in JHEP style - HYPER VERSION 



Quasiparticle Description of the QCD Plasma, 
Comparison with Lattice Results at Finite T and /i 



K. K. Szabo and A. I. Toth 

Institute for Theoretical Physics, Eotvos University, Pdzmdny 1, H- 11 17 Budapest, 
Hungary 



Abstract: We compare our 2 + 1 flavor, staggered QCD lattice results with a quasiparticle 
picture. We determine the pressure, the energy density, the baryon density, the speed of 
sound and the thermal masses as a function of T and /i b ■ For the available thermodynamic 
quantities the difference is a few percent between the results of the two approaches. We 
also give the phase diagram on the hb~T plane and estimate the critical chemical potential 
at vanishing temperature. 



Contents 

|l]. Introduction ^ 

|2|. Quasiparticle model ^ 

|3|. Comparison with lattice results ^ 

^ Conclusion 11 



1. Introduction 

QCD at high temperatures (T) and/or quark chemical potentials (/x) plays an important 
role in particle physics, since it describes relevant features of nature in extreme conditions. 
According to the standard picture of QCD, at high T and/or ^ there is a change from a 
state dominated by hadrons to a state dominated by partons (quarks and gluons). This 
transition happened in the early Universe (at essentially vanishing net baryon density Q) 
and probably happens in heavy ion collisions (at moderate but non-vanishing density) and 
in neutron stars (at large density, for which a rich phase structure is conjectured |2|, |3|, Q). 
There has been several studies in order to determine the thermodynamic properties of the 
QCD plasma. The weak coupling expansion of thermodynamic quantities shows a poor 
convergence near the transition temperature ||5|, |6|, |^] . Better convergence is expected from 
different improved techniques, e.g. screened perturbation theory |R |^, or a systematic 



rearrangement of the perturbative expansion |10|, or using an effective field theory in 3D 



1 11]. Note that noticeable deviations (~ 10%) from the ideal gas limit is expected upto 
temperatures as high as 1000 times the critical temperature (Tc). The lattice approach 
is applicable to investigate the plasma upto a several times Tc and at zero chemical po- 
tential |12, 13, 14|. The phenomenology of these lattice results has been widely analyzed 



in the literature. Successful quasiparticle descriptions [15, 16, 17| have been introduced 
to reproduce the properties of the QCD plasma. One can extrapolate these /U = re- 
sults to non-zero baryon densities by using the thermodynamical consistency of the model. 
Unfortunately, until recently it was not possible to compare the /i 7^ predictions of the 
quasiparticle model with direct lattice calculations. The lack of /i 7^ lattice results was a 
consequence of the sign problem. The fermionic determinant in the path integral becomes 
complex for non-vanishing //-s. This fact spoils any Monte-Carlo based technique used in 
numerical simulations. 
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Many suggestions were studied in detail to solve the sign problem and to give physical 
answers on the lattice at non-vanishing chemical potentials. Unfortunately, until recently 
none of them were successful. Recent interest in this field was initiated by the overlap 
improving multi-parameter reweighting method The phase diagram and the critical 

point were determined in the 2-1-1 flavor QCD on A'j = 4 lattices with staggered quarks 
1 19]. Several groups confirmed the results of Refs. |1^, 19 1 on the phase diagram 21, 22|. 



Furthermore, it became possible to determine the equation of state (EoS) at finite chemical 
potentials , too (for a recent review on lattice QCD at non- vanishing chemical potentials 
see Ref. [§4|). Thus, it would be of interest to see how these new lattice results can be 
described with the quasiparticle approach. This is the primary goal of the present paper. 

The paper is organized as follows. In Section 2 we briefly summarize the quasiparticle 
approach at = and show how to extend it to non-vanishing chemical potentials. Section 
3 compares our most recent lattice results with the prediction of the quasiparticle model. 
In Section 4 we summarize. 



2. Quasiparticle model 



In order to be self-contained we start with a brief review of the quasiparticle model ^ 
suggested in Ref. [16|. Consider a QCD plasma containing gluons {g), Ni number of 
"light" and number of "heavy" quarks (/ and h, respectively). These basic constituents 
have temperature and chemical potential dependent effective masses 



mj{T, fii) = + n*(T, fii), 



9,1, h, 



(2.1) 



where the moi-s are the rest masses, Hi-s are the chemical potentials of partons and the 
n*-s are the asymptotic values of the hard thermal/density-loop self-energies 



n* = 2ujq{mQ + u, 



1 



n* = - 

3 6 



Nr 



r2 + ^ 



(2.2) 
(2.3) 



where G^ is the effective gauge coupling (depending on T, /i, see below), and the summation 
is performed over the quark flavors (g). We are interested in the non-zero "light" quark 
chemical potential region therefore we set = for the "heavy" quarks. Note that 
the baryonic chemical potential is three times the "light" quark one. In our notations 
W = ^ = 

The thermodynamic potential the pressure p{T, /x) contains contribution from the ideal 
pressure of the quasiparticles (pl^) and from the pressure arising from their interactions 
(i?(r,M)): 



p{T,fi) = j;pf^(r,^,(^),mf) - B{T,is). 



(2.4) 



^When writing up this paper a different phenomenological approach was compared with lattice results 
at non- vanishing fi in Ref. [^. The quark-gluon plasma liquid model was found to be in close agreement 
with our lattice data. 
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The ideal gas pressure is the usual Fermi or Bose integral which takes into account the 
antiparticles in the quark pressure, too: 

(T, = A r de If ~ + - -n,), (2.5) 



iT,f^) = ^ de J'' . (2.6) 

Gvr^ ^m<,(T,A.) exp(e/r) - 1 

We do not fix the quark and gluon state multiplicities separately, but we demand that their 
ratio should be dg/{di^h) = 2 ■ 8/(2 ■ 3 ■ Ni^h). If we impose a stationarity condition |26| 



on the pressure then further thermodynamic quantities (such as energy density (e) , quark 
number density (nj) and entropy density (s)) will have only quasiparticle and mean field, 
B{T,fj,), contributions: 

e{T,fi)='£ei''[T,^,,{^^),mj]+B{T,^^), n,{T, fi) = nf^'iT, fii{fi),m^], (2.7) 
i 

where n^^ is the ideal quark number density. The functional form of the interaction 
pressure, B{T,fi), follows from the above condition, too. The derivatives of B{T,fi) in T 
and fi directions are easily accessible quantities 



dB 
df 



^ dmf dT' dfi ^~ ^ ^-^2 f,,, ' i 



dmf dfj, ' 



so B{T,fi) can be obtained by an appropriate line integral on the {T,iJ,) plane. 

What remains to be done is to determine the effective gauge coupling which appears 
in the formula of the self-energies. For vanishing chemical potentials it decreases logarith- 
mically with increasing temperature. A renormalization group inspired parametrization is 
as follows: 

G2(r,/i = 0) = ^ (2.9) 

[33-2(Ar,+Ar,)]iog(^) 

Imposing the Maxwell-relation between the derivatives of the quark number density and 
entropy 

dni 



ds 
5/i 



dT 



^ dm} dfi dmf dT ^ ' ' 



yields a first order, linear partial differential equation for G'^{T,fj,) with straightforwardly 
calculable aT,afj_,b coefficients 

ariT, fi, G^) ■ drG^ + a^{T, /i, G^) ■ d,,G^ = b{T, fi, G^). (2.11) 

This differential equation should be solved with the boundary condition at /i = (eq. 
( |2.9| )). Thus, the quasiparticle model is unambiguously defined by using thermodynamical 
consistency for non-zero /U-s above the critical line (Tc(/i)). Below Tc(//) the solution of 
the differential equation, thus the pressure is not unique. In this region the system has 
hadronic degrees of freedom instead of partonic. The quasiparticle model constructed from 
partons looses its validity. 
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3. Comparison with lattice results 



As we have already mentioned the sign problem of lattice QCD at finite chemical potentials 
spoils any Monte-Carlo method based on importance sampling. The recently proposed 
overlap ensuring multi-parameter reweighting method [|l8| enabled us to determine the 
EoS at non-zero temperatures and chemical potentials |^3|. The simulations were carried 
out on A't = 4 temporal and A^^ = 2A( . . . 3Nt spatial extension lattices with two "light" 
(rrii = 0.384rc) and one "heavy" (m/j ~ Tc) dynamical quarks. Note that the mass of 
the "heavy" quark corresponds approximately to the physical mass of the strange quark, 
whereas the mass of the "light" quarks is several times larger than the physical values for 
the up/down quarks. In the lattice analysis both the temperature and the baryon chemical 
potential covered the range upto /U^ ~ STc. 

Our lattice calculations were done on lattices with Nt = 4 temporal extension. In 
order to help the continuum interpretation we multiply the lattice results by the dominant 
T ^ oo correction factors between the Nt = 4 and the continuum case. We denote the 
T — > oo limit as the SB (Stefan-Boltzmann) case. 

_ P^^il^ = 0, T ^ oo, continuum) _ „ „ . ^ 

" p^^(^ = o,r^oo,A, = 4) - ^^-^^ 

^ _ Ap^^(^, T-.OC, continuum) _ ^^^^^ + O ( ^) . (3.2) 



Note that the m dependence of these factors are suppressed in the T — s- oo limit. Here 
Ap indicates the difference between the pressure at // / and the pressure at ^ = 0. The 
well-known continuum expressions for the non-interacting quark gluon plasma are 

p3^{^ = 0,T ^ oo, continuum) = [16 + 21(iV, + Nh)/2]ir'^ T"^ /90, (3.3) 
Ap^^{n,T^ oo, continuum) = Nii?T'^/2 + 0(/). (3.4) 

Including these dominant multiplicative corrections the results might be interpreted as 
continuum estimates. (Clearly, the appropriate - but more CPU-consuming - way is to 
carry out the lattice calculations at A^ > 4 and perform an extrapolation to A^ — > oo.) 

The quasiparticle model has four free parameters: A and to determine the gauge 
coupling, dg the gluon multiplicity and B{Tc) the integration constant in the interaction 
pressure. They were used as fit parameters in order to receive the least possible difference 
between the thermodynamic observables measured on the lattice and the ones predicted 
by the quasiparticle description. In principle it is sufficient to use only the lattice results 
at vanishing in the fitting procedure since afterwards the thermodynamical consistency 
unambigously defines the quasiparticle approach; but this could lead to large differences 
between the two types of descriptions at non-vanishing //-s. It turned out to be more 
sensible to use both the fj, = and the /U 7^ lattice data for fitting in order to gain 
a better agreement between the two approaches for higher ^-s. Therefore we fitted the 
quasiparticle picture on the (e — 3p)/T^ at = and Ap/T'^ at our highest chemical 
potential {fiB ~ 490 MeV) lattice results, simultaneously. The chi-square function (x^) was 
constructed by using the statistical errors of the lattice data. We constructed statistical 
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errors for the fit parameters, which are to give the 68.3% confidence level corresponding to 
the bounds of the Ax^ = 4.72 interval |27|. Clearly, the parameter region according to this 



confidence level has a complicated shape in the four-dimensional parameter space. The 
quoted errors just indicate the borders of this domain. Furthermore, these errors should 
be taken with a grain of salt, since the constructed functions contains no information 
about the systematical uncertainties of our Nt = 4: lattice data. We found it to be more 
reliable to fit on (e — 3p)/T^ lattice results instead of p/T'^ at /i = 0, since minimizing the 
function of the pressure yielded an interaction measure without a turning-point around 
1.3Tc, which is typical of the lattice results. The best parameter values with statistical 
uncertainties can be found in Table |l[ 



A 


10.1 


Ts/T, 


-0.85 


dg 


16.4 




174.3 til MeV 



Table 1: The best fit parameters of the quasiparticle model to our lattice results. 



Using the best fit parameters we can obtain the temperature and chemical potential 
dependence of several thermodynamic quantities and compare them with lattice data. Note 
that we use Tc = 170 MeV as the overall scale in the results. In formulas and figures Tc 
is used as the transition temperature at vanishing chemical potential, whereas Tc{^b) 
indicates the transition temperature at non-vanishing chemical potentials. We chose the 
chemical potential values so as to cover the characteristic /i-s occuring in the heavy ion 
collisions. At our smallest chemical potential value, 27 MeV, and at around 300 MeV 
Au-t-Au colliding experiments were carried out at RHIC [p8| , p9| , whereas at 253 MeV 
Pb-hPb collisions took place at SPS |0|. 

Figures ^ and |2| show the behavior of the "interaction measure", (e — 3p), and the 
pressure (both normalized by T^) at vanishing chemical potentials. Both the lattice result 
and the prediction of the quasiparticle approach are given. On the one hand there is a good 
agreement between the two techniques for (e — 3p), on the other hand for the pressure one 
finds a larger difference. It is easy to understand the reason for that. First of all we used 
(e — 3p) in our fitting procedure, therefore a better agreement is expected for this quantity 
than for p, which is a predicition of the model. Secondly, the pressure can be obtained as 
an integral of the interaction measure 



rp4 



dT e — 3p 

rp rp4 



(3.5) 



Differences in (e — 3p) at small temperatures dominates the above integral. Though the 
difference between the lattice and the quasiparticle results for (e — 3p) is rather small at 
large T one observes a 20% difference around Tc. Therefore we end up with an ~ 10% 
difference for p even in the large T region. There might be several explanations for the 
differences between the two approaches. It can be that the lattice results in the continuum 
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Figure 1: The interaction measure, (e — 3p), normalized by T'* as a function of T/Tc at /i = 0. 
Tire line corresponds to the quasiparticle model, the points are lattice results multiplied by Cp. The 
errorbars show the statistical uncertainties. 

T[MeV] 
200 300 400 500 




T/T„ 

Figure 2: The pressure normalized by as a function of T/Tc at /i = 0. The line corresponds 
to the quasiparticle model, the points are lattice results after a multiplication by Cp. The arrow 
indicates the high temperature ideal gas pressure of the QCD plasma (SB Hmit). 

limit {Nt oo) change the picture and lead to smaller discrepancies. Alternatively, it 
also can be that the quasiparticle approach simplifies the interaction in the quark-gluon 
plasma. This means that the predicition of this method can be ~ 10% off. 

We show the temperature dependence of Ap/T^ in Figure ^ and that of ub/T^ in 
Figure |. Both fi gures indicate a nice agreement between the two types of QCD plasma 
descriptions (the largest difference is ~ 5% at low temperatures and high chemical poten- 
tial). The temperature dependence of the baryon density is similar to that of the pressure 
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Figure 3: Ap = p{fi ^ 0, T) -p{fi = 0, T) normalized by T'^ as a function of T/T^ for fig = 27, 45 
and 60 MeV in the left panel and for fis = 140,290,410 and 490 MeV in the right panel (upper 
curves correspond to larger chemical potentials). The lines represent the quasiparticle model, the 
points are lattice results multiplied by c^. The statistical uncertainties are smaller then the point 
size. 



T[MeV] TLMevj 
200 300 400 500 200 300 400 500 




Figure 4: The baryon number density normalized by T'^ as a function of T/Tc for jiB ^ 27,45 
and 60 MeV in the left panel and for = 140,290,410 and 490 MeV in the right panel (upper 
curves correspond to larger chemical potentials). The lines represent the quasiparticle model, the 
points are lattice results multiplied by c^. The statistical uncertainties are smaller then the point 
size. 



since they are connected by the formula 



d 



T 



riB 
T3- 



(3.6) 
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Figure 5: Ap of the interacting QCD plasma normalized by Ap^^ of the free gas (SB) as a 
function of T/Tc for fiB = 60,140,290,410 and 490 MeV. The lines represent the quasiparticle 
model, the points are lattice results. The statistical uncertainties are smaller then the point size. 

We observe an interesting scaling behavior of Ap normalized by Ap^^ (Figure |5|). This 
quantity {Ap/Ap^^) depends only on T/Tc and it is almost independent of the chemical 
potential. This scaling behaviour is less accurate around T/Tc and it gets more and more 
precise for higher temperatures. 

It is of extreme importance to determine the phase line (the line, which separates the 
phases dominated by hadrons and partons). One is particularly interested in results at low 
temperatures, for which only model estimates are available. First we use the quasiparticle 
picture, then an extrapolation to the lattice data in order to determine the phase line. 

We defined a quasiparticle critical line based on the following observation (see Figure|2|) : 
the pressure around Tc is a lot lower than the pressure at high temperatures, which is in 
the order of p^^. A simple explanation of this fact is that the number of degrees of 
freedom in the QGP is much bigger than in the hadronic phase, which is a dilute gas of 
pions. So we set p to be a few percent of p^^ along the phase line (note that the result 
is rather insensitive to the exact value of this percentage as illustrated in Figure P). The 
line defined above is depicted also at rather low temperatures and we get ~ 1800 MeV for 
the point of intersection on the T = axes. We must interpret this value as the critical 
chemical potential {^%^^) with great care since for the low temperature, high chemical 
potential region the color-superconducting phase is conjectured [^, which is obviously 
not contained in this quasiparticle framework. The quasiparticle transition line shows a 
nice agreement with the directly measured lattice transition line in the hb < "^Tc region 
(Figure |6|). Outside of this region we do not have direct lattice results. It is intriguing to 
see the behavior of the critical lines when T ^ 0. We extended our lattice transition line 
to zero temperature by using two fits namely a second and a fourth order polynom. Due to 
versus symmetry we keep only even order terms in the fitting procedure. The bands 
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Figure 6: The critical lines on the temperature versus chemical potential plane. The solid curves 
are quasiparticle critical lines which were obtained with the best fit parameters. The lower curve 
corresponds to the pressure equals to zero criteria while the pressure equals to 0.14^"^^ along the 
upper curve. The dashed and dashed dotted (lower and upper) lines are the extrapolations of 
the lattice critical line to T = using a 4th and a 2nd order polynom, respectively. The dotted 
lines show the influence of the statistical uncertainties of the lattice results. Note that our critical 
lines loose their validity towards higher chemical potentials, where the color superconducting phase 
comes into play. 
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Figure 7: The speed of sound squared in the QCD plasma as a function of T/Tc{^b) at various 
chemical potentials (/is = 0, 900, 1200 and 1500 MeV). Upper curves correspond to larger values of 
/is. The arrow indicates the speed of sound squared in the ideal gas limit ((cf^)^ ~ V^)- 

in Figure ^ indicate the statistical uncertainties of the fits. The two curves agree nicely 
in the directly measured fiB < 3Tc region. They deviate for smaller temperatures and 
larger chemical potentials. The quadratic polynomial fit predicts at vanishing temperature 
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Figure 8: The coupling constant squared as a function of T/Tc for various chemical potentials 
{hb = 0,900, 1200 and 1500 MeV). Upper curves correspond to smaller values of ^b- 

/ig'* 2000 MeV (which exceeds expectations). The quartic approximation gives for the 
same quantity /i™* ~ 1400 MeV. The difference between the predicitions of the two kinds 
of polynoms might be interpreted as a rough estimate for the systematic uncertainty of the 
extrapolation to low temperatures. 

The dynamical properties of the plasma phase is primarily determined by the speed of 
sound 

c?(r,/.B) = |. (3.7) 

As it can be seen in Figure |^ increasing chemical potential yields higher speed of sound. 

One can also look for the predictions of the quasiparticle model for quantities which 
are straightforward predicitions of the quasiparticle approach however, it is difficult to 
obtain them directly on the lattice. Such important quantity is the gauge coupling of the 
quasiparticle picture, which is shown in Figure ^. Note that the static potential at finite 
chemical potentials can be in principle measured on the lattice. Using the lattice potential 
one can easily define the gauge coupling. Clearly, the parametrization of the gauge coupling 
in the quasiparticle picture is not equivalent with the above mentioned coupling. 

In order to understand the typical degrees of freedom in the high temperature QCD 
phase at non-vanishing chemical potentials it is instructive to study the effective masses of 
the quasiparticles. Figure ^ shows the temperature dependence of the light quark masses for 
different chemical potentials, whereas Figure |l^ presents the results on the gluon masses. 
According to our observations these quantities are almost independent of the chemical 
potential in the fJ-B < 3Tc region. Due to the stationarity condition dp/dm\ = small 
changes in the mass do not change the pressure. It means that the chemical potential 
dependence of the thermodynamic quantities are primarily coming from the direct 
dependence of the Fermi integrals. 
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Figure 9: The effective mass of the light quarks as a function of T/Tc for various chemical 
potentials (a^s = 0, 900, 1200 and 1500 MeV, right to left). 
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Figure 10: The effective mass of the gluon as a function of T/Tc for various chemical potentials 
{fjLB = 0,900, 1200 and 1500 MeV). Upper curves correspond to smaller values of ^b- 

4. Conclusion 

In this paper we studied the quasiparticle approach to describe the equation of state of 
the hot QCD plasma. We have fitted the free parameters of the model by using our 2 + 1 
flavor, dynamical staggered QCD lattice results in the region T,hb < STc. After calculat- 
ing the pressure, interaction measure and density we found a good agreement between the 
quasiparticle predictions and our lattice data. The model successfully justifies the scaling 
behavior of Ap/Ap^^ observed in the lattice calculations. We gave some confidence in- 
tervals for the fit parameters in spite of the lack of continuum extrapolations. Using the 
best fit parameters the quasiparticle critical line and the speed of sound were given for 
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higher values of hb- The zero temperature hmit of the quasiparticle and lattice critical 
lines (critical /ib value) cover the //g^* « 1400. . . 2000 MeV region. 
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